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' Abstract 

The generalization of the Jaynes-Cummings (GJC) Model is proposed. In this 
model, the electromagnetic radiation is described by a Hamiltonian generalizing the 
harmonic oscillator to take into account some nonlinear effects which can occurs in 
the experimental situations. The dynamical superalgebra and supercoherent states 
of the related model are explicitly constructed. A relevant quantities (total number 
\Q [ of particles, energy and atomic inversion) are computed. 

o 
o 
o 

1 Introduction 

The Jaynes-Cummings (JC) model [1] which is extensively used in quantum optics 



describes, in its simplest version, the interaction of a cavity mode with two-level system 

Hj C = u(a + a~ + ^)cr + yff 3 + n(a + a_ + a~a + ), (1) 

where a + and a~ are the photon creation and annihilation operators, a± = |(o~i ± 02), 
with o"i, 02 and 03 are the Pauli matrices and <tq is the identity matrix. Moreover, k is 
a coupling constant, u is the radiatif field mode frequency and lj the atomic frequency. 
The interest of this model, its solvability and its applications, has long been discussed 
[1, 2, 3]. Over the last two decades, there has been intensive study [4 and references 
quoted therein] on the solvable Jaynes-Cummings model and its various extensions, such 
as intensity depending coupling constants, two photons or multiphoton transitions and 
two or three cavity modes for three-level atoms. These models have found their applica- 
tions in laser trapping and cooling atoms [5] and quantum nondemolition measurements 
[6]. Furthermore, the Jaynes-Cummings model constitutes now the basis for a vast ar- 
ray of the current experiments on fundations of quantum mechanics involved entangled 
states [7, and references therein]. In other hand, the supergroup theoretical approach 
to Jaynes-Cummings model has opened the way to relate the exact solvability of this 
model and representation theory of super algebras. Indeed the Hamiltonian Hjc is an 
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element of the u(l/l) superalgebra [8]. In the absence of coupling (k = 0) and for ex- 
act resonance (cu = cu ), the u(l/l) dynamical superalgebra reduces to si (1/1) one and 
the JC model coincides with the supersymmetric harmonic oscillator. More recently, the 
investigations of a class of shape-invariant bound state problem, which represents two- 
level system, leads to generalized Jaynes-Cummings model [9, 10, 11]. In the case of the 
simplest shape-invariant system, namely the harmonic oscillator, the generalized Jaynes- 
Cummings model reduces to standard one. 

In this paper we shall address the generalization and quantum characteristics of the 
Jaynes-Cummings model. Besides the eigenvalues and eigenvectors, we give the super- 
coherent states of the related model. It is found that the generalized Jaynes-Cummings 
model is governed by a nonlinear superalgebra u(l/l) which reduces to well known u(l/l) 
occurring in the standard (JC) model [12, 13]. We compute the total number of photons 
and the energy. We find that the atomic inversion exhibits Rabi oscillations. 

The paper is organized as follows: In Section 2, we introduce the generalized super- 
symmetric quantum oscillator and we construct the corresponding supercoherent states. 
Exact spectrum of the generalized Jaynes-Cummings model is given in section 3. Section 
4 is devoted to nonlinear dynamical superalgebra u(l/l) of the (GJC) model which is 
useful to construct the supercoherent states adapted to our model (Section 5). Using the 
latter set of super-states, we compute in section 6 some relevant physical quantities. The 
last section concerns the conclusion of this work. 



2 Generalized Supersymmetric Quantum Oscillators 

We begin by introducing the generalized supersymmetric quantum oscillators. Let us 
Consider a Hamiltonian H with a discrete spectrum which is bound below and has been 
adjusted so that H > 0. We assume that the eigenstates of H are non-degenerate. The 
eigenstates |\I/ n ) of H are orthonormal vectors and they satisfy 



In a general setting, we also assume that the energies e , e±, e 2 , ... are positive and verify 
e n+i > e n- The ground state energy is e = 0. We define the creation and the annihilation 
operators A + and A~ , respectively, such that the Hamiltonian can be factorized as 



H\V n ) = e n |* n ). 



(2) 



H = A + A~. 



(3) 



The action of the operators A + and A 



on the states \^ n ) are given by 



A + \^ n ) 
A~\* n ) 



(e n+1 )2|^ n+1 ) 
(e n )5 



(4) 



implemented by the action of A on the gound state |^ ) 



A-|* ) = 0. 



(5) 



The commutator of A + and A is defined by 



[A-,A + ]=G(N), (6) 
where the operator G(N) is defined through his action on \^/ n ) 

G(N)\^ n ) = (e n+1 -e n )\^ n ). (7) 
We define the number operator N as 

N\V n )=n\V n ), (8) 



(9) 



N is in general different from the product A + A (=H). We can see that the number 
operator satisfy 

[A+,N] = A+ 
[A~,N] = -A-. 

Here, we consider two generalized oscillators systems which has extensively studied in the 
literature. The first concerns the so-called generalized deformed oscillator [14] and the 
second one is the x 4 -anharmonic oscillator [15]. The physical interests of this two systems 
have been extensively enumerated [see the references quoted in [14, 15]]. Here, we recall 
their eigenstates and eigenvalues to construct the supersymmetric generalized quantum 
oscillator and the corresponding supercoherent states. 



To introduce the generalized deformed oscillator, the procedure of [14] requires the 
existence of a map from the usual harmonic oscillator algebra generated by annihilation 
and creation operators a~and a + satisfying the standard canonical commutation relations, 
to the new one generated by A~ and A + 

A- = a-f(N) A + = f(N)a + (10) 

N being the number operator N = a + a~ and the function / is given by 

f(N) = N + m (11) 

The Hamiltonian of the obtained generalized harmonic oscillator is then given by 

H = A + A~ = N(N + m) (12) 

with eigenvalues 

e n = n(n + m) (13) 

The Fock states |\l/ n ) = \n, m) are labelled by the integers m and n — 0, 1, 2, .... 
It is clear that the operator G(N) (7), in this case is given by 

G(N) = 2N + m+l. (14) 

Then, the operators A + ,A~ and G(N) satisfy the relations (6) and (9). Note that other 
choices of the function / are possible. We remark also that when f(N) = 1, we have the 
ordinary harmonic oscillator. 



The other nonlinear oscillator that we consider is the x 4 -anharmonic oscillator. The 
Hamiltonian, describing this system, is 



H = a + a~ + |(a" + a+) 4 - 5 (15) 

where a + and a~ are the creation and annihilation operators for the harmonic oscillator. 
The quantity 5 is given by 

H £ -f e2 (16) 

which vanishes when e = and H reduces in this limit to the standard harmonic oscillator 
Hamiltonian. The Hamiltonian H can be factorized in the following form [15] 

H = A + A- (17) 

in terms of A + and A~ which are expressed as some functions of a + and a~ (for the 
expressions of these functions see[15]). The energy levels are given by 

e n = n + |e(n 2 + n), n = 0,1,2,... (18) 

The positive parameter e can be seen as one taking into account some non linearity of the 
radiatif field arising from some perturbative effects occurring in experimental situations. 
Note that we keep terms only up to e which is the standard first-order perturbation result. 
The Hilbert space of this system is easily constructed in the same way as the standard 
harmonic oscillator. It is spanned by the states \n, e), n — 0, 1, 2, ... which is generated by 
the action of A + on the ground state |0, e). The operator G(N) is 

G(N) = 1 + 3e(n+ 1) (19) 

Here, again one can verify that the relations (6) and (9) are satisfied by the creation and 
annihilation operators corresponding to the x 4 -anharmonic system. 

In supersymmetric quantum mechanics, one consider the so-called supersymmetric 
Hamiltonian which is defined by 

H ( A ~ A+ \-( H + M m\ 



A+A- J V H 

where H+ = A~A + and H = A + A~ = H are the so-called supersymmetric partner 
Hamiltonians. 

Working in the Hilbert space 

h = h (g) h f = -) = ( , ° ) , |* n , +) = ( l ^ ) ; n = 0, 1, 2, ...} , (21) 



the eigenstates are 
I0o> 




l*o) 



(22) 



with the energies E = and E n>0 = e n . Because, we are interested by generalized 
quantum oscillator, the states |^ n ) are \n,m) for the generalized deformed oscillator and 



\^ n ) are \n,e) for the x 4 -anharmonic oscillator. As we have mentioned previously these 
two quantum systems will be used to extend the JC model and we will compute some 
relevant physical quantities, like the mean values of the total number operators, the energy 
and the atomic inversion, over the coherent states of the (GJC) model. The latter will 
be obtained from the supercoherent states corresponding the generalized supersymmetric 
oscillator. So for H susy , we consider the supercoherent states (linear combination of the 
fondamental coherent states) 

|.,/5) = cos^,M+sin^f^V (23) 



2 \ \z)J 2 V 

with j3 = |e i$ and \z) are the coherent states, corresponding to H_, defined by [16, 17, 
18, 19] 

oo n 

I*> = K(M)E7TT7tI*»>> (24) 
n =o (e(n))a 

where e(n) = e\...e n for n = 1,2,... and we set e(0) = 1. The normalization constant 
N(|z|) is calculated from the normalization condition (z\z) = 1 and it is given by 

oo I 1 2n 

(%I))- 2 = Ett- (25) 

e ( n ) 

Let us mention that the coherent states for x 4 -anharmonic oscillator has been studied in 
[19]. They are given by 

oo nn 

\z) = a(\z\)Y (——-, — 2r)^ n |n,e), (26) 

with 

« = -ffi±^. (27) 

For the deformed generalized harmonic oscillator, we construct the coherent states in 
the same way that one which gives Barut-Girardello coherent states of the su(l, 1) algebra 
[20]. Note that the algebra generated by {A + , A~,G(N) = 2N + m + 1} is isomorphic to 
su(l, 1) ~ s/(2, (R)) ~ so(2, 1). Indeed, the creation A + and annihilation A~ operators 
satisfy the following commutation relations 

[A-,A+] = G(N), [A±,G(N)] = T A±. (28) 

A more familiar basis for su(l, 1) algebra is given by 

■J- = ^A-, J+ = ^A+, J 12 = 1G(N), (29) 

with the following commutation relations 

[J-,J+] = Ji2, [J±,M = tJ±- (30) 

Barut and Girardello introduced the su(l, 1) coherent states as eigenvectors of Jl [20] 

J-\z) = z\z) 

\z) = b(\z\) g {Vlz)n r \n,m). (31) 

n=0 (n!r(n+m+l))2 



The normalization constant is given by 



m = (r ( m + i))i 

Where the oFi(m + l,2|z| 2 ) is the hypergeometric function. The coherent states (26) ( 
for the a; 4 -anharmonic oscillator) and (31) (fot the generalized deformed oscillator) can 
be obtained simply from equation (24) by replacing in the expressions of e(n) the energies 
by their corresponding values for each considered system. Note that the coherent states 
(26) and (31) will be useful to build up ones of (GJC) model (see section 5). Remark also 
that the resolution to identity of such states has not been discussed here. However, the 
measures by respect with the previous sets of coherent states are over-completes can be 
computed in a very easy way following the approach developed in [13, 17]. 



3 Eigenstates and Eigenvalues 

To generalize the JC model, let us consider the Hamiltonian 

H GJC = A + } + -u [A~, A + ]a 3 + n(A + a^ + A~a+). (33) 

This last expression generalizes the ordinary JC Hamiltonian. In fact, when the creation 
and the annihilation operators are those associated with the harmonic oscillator, the 
above Hamiltonian reduces to the well known JC model (Eq.(l)). We note also that the 
Hamiltonian Hgjc is supersymmetric when uj = ojq (exact resonance) and k — (absence 
of coupling) 

H gjc (k = 0,u = u ) = {Q-,Q + }, (34) 
where the supercharges operators are defined by 

Q~ = iy/uia + <j-, Q + = —iyfujar cr + , (35) 

and satisfy the relations 

(Q-) 2 = 0, (Q+) 2 = 0, [0*^ = 0. (36) 

We note that the generalized Hamiltonian of JC model (33) is different from ones con- 
sidered in the references [21, 22, 23]. The main purpose of this section is to show that 
generalized JC model can be solved analytically. The solutions become ones correspond- 
ing to JC model when A + = a + and A" = a~ (i.e. harmonic oscillator). 



The diagonalization of the Hamiltonian Hgjc is easily carried out (in the same way 
that the standard JC model) and leads to the eigenstates 

\Eo) = |*o, -> 

\ E n) = T ^ W) [S(n + l)\y n ,+)-Q(n+l)\y n+1 ,-}] (37) 
\E~ +1 ) = p^[Q(n+l)|tt n ,+) + S(n + l)|tt n+1 ,-)], 



where 

Q(n + 1) = «(e n+ i)2 

S(n + 1) = ^(e n+1 + ^(^^) 2 )^ + ^^|^, (38) 

P(n) = ((S(n)) 2 + (Q(n)) 2 )* 

with A_ = uj — ujq. 



The corresponding eigenvalues are given by 



E; 1+1 = ^ eri+1 + ^(e n+2 + e n ) + (« 2 e n+1 + %(^|^) 2 )i 



(39) 



where A + = uj + u . 

Interesting particular cases arise from the former results. We start with the exact 
resonance (A_ = 0). In this case, the Hamiltonian H GJC takes the form 

H GJC (A_ = 0) = \oj{A~, A + } + A + }a 3 + k{A + o_ + A~a+). (40) 

The eigenstates of the resonant generalized JC model are then 

\K) = ^[|*n,+>-|*n+l,->] ur , 
\E~ + l) = ^[|*n,+> + |*n + l,->]- 

and the corresponding eigenvalues are 

E+ = ue n+1 - /t v /e^T , 42 , 
E n+i = ue n+1 + Ky/e n+1 . 

The resonant generalized JC model is reduced to supersymmetric Hamiltonian when the 
coupling constant k vanishes 

H susy = H GJC (A_ = = «) = u ( J A _ ) . (43) 

Finally, we remark that when e n = n (spectrum of the quantized radiatif field), the 
eigenstates and eigenvalues (37) and (39) coincide with ones corresponding to standard 
JC model [12] (see also [13]). 

4 Dynamical Superalgebra 

Like standard JC model [12], the generalized JC model can be written as a linear 
combination of two even operators Ni and N 2 and two odd operators Q~ and Q + ; 

H G jc = ^N, - ^-N 2 + -£=(Q+ - Q~), (44) 

Z Z \ /UJ 



(45) 



where 

A~A+ \ „ f-A+A- 



lVj = 1 .1.1 J" Xl { -A-A+ 
Q~ = i^fujA+O-, Q + = —i^JujA~a + . 



These operators satisfy the following relations 

[N u N 2 ] = 0, [N l7 Q~\ = [N ± , Q + \ = 0, 

[N 2 , Q + ] = {Q + 9n + 9nQ + ), [N 2 , Q-] = -{Q-g N + 9nQ-), (46) 
{Q-,Q+} = uN 1 , (Q-y = (Q+) 2 = 0, 

where the even operator g N is defined by 

g N = G(N)a . (47) 

The algebra generated by {N lt N 2 , Q~ , Q + } can be seen as a non linear (or deformed) 
version of the superalgebra u(l/l). Indeed, when G(N) = 1 (a situation which occurs in 
the ordinary JC model), we have the w(l/l) superalgebra. In the situation, in which we 
are in presence of an exact resonance w = w , the generalized JC model is written as 

H GJC = cuN, + ^=(Q+-Q-), (48) 

in terms of the operators N±, Q + and Q~ which satisfy the structural relations of the 
superalgebra s/(l/l). It is also easy to see that for uj = u and k = 0, we have 

H gjc = ujN 1 = {Q+,Q-}, 

[Q ± ,H ajc ] = o, (Q-) 2 = (g + ) 2 = o. 

5 Supercoherent states for Hqjc 

To construct the supercoherent states for Hgjc, we start by looking for an unitary 
operator U which connects the Hamiltonians Hp (diagonal in the fondamental space of 
states {|* n , +), |* n , -)}) and H GJC 

H D = U+H GJC U=( H + * ), (50) 

where 



H + = ^g(N + 1) + ^(g(N + 2) + g(N)) - (^g(N + 1) + ^n + 2)- 9 (n) )2) i 
H- = ^g(N) + ^t(g(N + l)+g(N-l)) + (^g(N) + ^( rfW-sW )^. 

(51) 

The operators g(N + k) (k — — 1, 0, 1, 2) are defined by 

g(N + k)\* n ) = e n+fc |^ n ). (52) 

The operator U takes the following form 

u= ^S(N + l) r^A-\ (53) 
\ 71 p(n+i) p(N)°y iy i / 



where 



S(N) = ± « N+1 >- « N - » + .fa (N) + ^ N+1) - - lj ))*, (54) 



and 

P(N) = (S(N) 2 + K 2 g(N))*. (55) 
The operator U can be written as follows 

U = e- z , (56) 

with 

Z = A + h(N+ l)tr_ - h(N+ l)A~a+. (57) 
Developing e~ z and using the following identities of the Hermitian skew operator Z, 

Z 2 ? = (-mh(N + l)) 2p (g(N + l))V + (7_ + (-iy(h(N)) 2 *>(g(N)ya-a + (58) 

and 

Z 2 P+ i = (-iy+\h(N + l)) 2p+1 (g(N + l)) p A-a + + (-l) p A + (h(N + l)) 2p+1 (g(N + 1))V_. 

(59) 

It becomes that [/ is given 



/ cos(/i(JV + l)(^(iV + l))3) — ^^ sin (/ l (AT + l) 1 — T )A- 



(60) 

where /i(iV) = 7 l arctan(/t^S^P). Of course, for g(N) = N (standard JC model), we 

y/g(N) 

recover the results obtained in [12]. 

Using the unitary operator U, we introduce the coherent states for generalized JC 
Hamiltonian as follows 

\z,p) GJC = U\z,P) 

= ^(N)E^=(sin|e^) + cosfK-)). ( 61 ) 

n=0 V ' 

To compute some relevant physical quantities of the generalized JC model, we consider 
the time evolution of states (61) 



\z,/3,t) G jc = e- itH <"°\z,/3). (62) 



We get 



oo n n a 

\z,P,t)QJC = *(M) E 7T^( sin 9 el% ~ UEn \ E n)+^s-e- itE ^\E-)). (63) 

Using the latter states, we will compute the average value of the total value of photons, 
energy and atomic inversion. 



6 Total number of photons, energy and atomic inver- 
sion 



6.1 Total number of particules 

The mean values of the operator N over the states (63) is given by 

oo I _ I In a n 

(AO = K(M) 2 E Vr^ 2 |«I*K> + ™s 2 -(E~\N\E-)). (64) 

n=0 e \ n ) Z Z 

A direct computation of matrix elements occuring in the last expression gives 

oo I |2n n n 

( N ) = «(l^) 2 E l zL(^ 2 ^n + i + cos 2 -en). (65) 

n=0 e \ n ) Z Z 

where the energies e n are given (13) (resp.(18)) for the generalized deformed oscillator 
(resp. for the x 4 -anharmonic system). 

6.2 Energy 

The energy in coherent states (63) is, simply, given by 

oo I 1 2n n n 

(H G JC) = K(M) 2 E VT( Sin2 + C ° S2 E n)' ( 66 ) 

n=0 e l"/ Z Z 

The i? n + and E n ~ are given by Eqs (39). 

6.3 Atomic inversion 

First, if we start with supercoherent states (61), we see that the average value of the 
third component of the spin is time independent. However, as it is well known, if the radi- 
atif field is prepared in a coherent state, the atomic inversion consists of Rabi oscillations. 
The temporal dependence of (<7 3 ) appears when we use the generalized supercoherent 
states (63). Indeed, we get 



<^3> = W++ + (<7 3 )— + (<7 3 )+- + (<7 3 )-+, (67) 



with 



<^3> ++ = i(l-cosWN) 2 £ ^q(n+l) 



n=0 



<"3>- = -|(l + cos^)N(|^|) 2 g^gH 
(<7 3 )+_ = isin^(l^l) 2 £ ^i e i '( 2KS ( n+1 )+*- Q ) 



(68) 



where 



2 D± " n ^ Q e(n) " s(n+l) 

(a 3 )_ + = Isin^(|,|) 2 E^e-^(^H^)^, 

n=0 v ' v ' 

g(» + l)= / ^ e _" +2 - e " } , q(0) = l 

V e »+l + ( — ) 2 ( e «+2- e ») 2 

s(n + l) = ^ j-i , s(0) = l (69) 

V e « + l + (-4^) 2 ( e " + 2- e n) 2 

ct = arctanz, 



The result of the computation of (03) shows that the time dependance comes from the 
value (0-3)+- + (03)-+ and we obtain 

. |^| 2 " +1 .sin^cos(t(2/ts(n+ 1) + $-«)). 

n=o e ( n ) s(n + l) 

This expression have an oscillating behaviour that characterizes the atomic inversion in 
the generalized JC model. 

7 Conclusion 

In this work, we developed the generalization of the Jaynes-Cummings model where 
the radiatif field is replaced by generalized harmonic oscillators. We shown the exact solv- 
ability of the model. The role of the nonlinear dynamical superalgebra u(l/l), governing 
the evolution of the related model, is important in the construction of the supercoherent 
states over which we compute the energy, average number of photons and the atomic 
inversion. These states has been constructed using the unitary transformations expressed 
in terms of the generators of the u(l/l) superalgebra. Furthermore, it is shown that the 
time dependent supercoherent states have the advantage to obtain the Rabi oscillations. 

Finally, we note that the importance of the generalized Jaynes-Cummings model is 
due not only to its exact solvability but also arises from its quantum effects such as revival 
of atomic inversion. In our opinion the generalization and results presented here can be 
used to give a realistic description of the nonlinear process of the interaction of an atom 
and radiation field. Clearly, only a confrontation with experimental measures can valid 
the results of this work. 
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